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1  echnical  Report  Summary 

The  research  undertaken  in  this  project  is  directed  to  three 
principal  tasks.  The  first  of  these  is  the  application  of  three-dimensional 
seismic  ray  tracing  to  problems  of  hypocenter  location  and  focal  mechanism 
dete  ruination  thereby  freeing  the  analyst  from  the  unrealistic  assumptions 
of  a  spherically  symmetric  earth.  One  of  the  products  of  this  task  will  be 
a  greatly  improved  model  of  upper  mantle  structure  and  wave  propagation 
under  the  Hindu  Kush  region  of  central  Asia.  The  second  major  task  is 
generally  directed  to  exploring  the  full  capabilities  offered  by  the  High 
Gain  Long  Period  network  for  discrimination  of  earthquakes  and  explosions, 
especially  at  low  energy  levels.  The  third  task  as  an  attempt  to  establish 
a  classification  of  Asian  earthquakes  in  interesting  source  regions  on  the 

basis  of  spectral  content,  dominant  focal  mechanism,  b-  values,  and  other 
descriptors. 

The  principal  achievements  to  date  under  the  first  task  are  the 
completion  of  the  development  of  a  computer  code  to  treat  ray  tracing 
in  a  generally  heterogeneous  medium  and  the  successful  testing  of  this 
code  by  application  to  the  central  part  of  the  Aleutian  arc.  Significant 
improvement  in  hypocenter  locations  and  greatly  increased  ability  to 
relate  the  events  to  the  major  structural  elements  have  been  achieved. 

The  effects  of  upper  mantle  structure  on  focal  mechanism  solutions  and 
dT/dA  measurements  at  teJeseismic  distances  has  also  been  studied, 
with  results  presented  in  this  report. 
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An  at  related  effort  under  th'.s  task  is  the  investigation  of  mean 
P  travel-time  residuals  at  stations  of  the  Soviet  seismographic  network. 
These  results  will  be  used  as  part  of  the  input  for  the  development  of 
models  of  crustal  and  upper  mantle  structures  in  central  Asia. 

Work  on  long -period  body  waves  has  developed  along  two  lines. 
Spectra  of  body  waves  from  c  ustal  events  recorded  on  standard  long- 
period  instruments  of  the  \VW  SSN  have  been  interpreted  in  terms  of  focal 
depth.  Data  from  the  HCLP  network  will  be  used  for  this  purpose  also, 
but  that  system  is  not  yet  up  to  full  operational  status.  Studies  to  date 
of  explosion  and  earthquake  seismograms  from  the  HGLP  netwoi  .c  have 
revealed  that  long-period  P  waves  are  essentailly  absent  of  explosions, 
but  SV  type  body  waves  are  strongly  generated  and  are  the  only  body  wave 
type  clearly  seen  on  the  HGLP  records. 

Under  the  thud  task,  P  and  S  wave  spectra,  magnitudes  (long  and 
short-period  body  wave  magnitudes  and  surface  wave  magnitudes),  and 
focal  mechanisms  have  been  derived  for  events  occuring  in  from  selected 
source  regions  of  central  Asia.  This  task  is  still  in  the  data -gathering 
and  basic  analysis  stage  and  no  general  conclusions  have  been  reached. 
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jj _ Applications  of  Seismic  Kny  f  rue  ins: 

E.R.  Fngdal.l  ' 

Tlic  problem  >f  ray  tracing  in  a  generally  heterogeneous  medium 
has  been  treated  using  the  calculus  of  variations  and  Fermat's  principle 
of  stationary  time.  The  solution  in  geocentric  spheric?)  polar  coordi¬ 
nates  is  a  system  of  five  simultaneous  first-order  differential  equations 
giving  the  variation  with  time  of  the  position  (r,  9,  0)  of  a  point  on  a  ray 
and  the  direction  of  propagation  (  i  =  incident  angle,  a  =  azimuth)  from 

this  point  in  terms  of  the  wave  speed  (v)  and  its  spatial  derivatives  in 
the  medium  , 

A  computer  program  has  been  developed  to  treat  this  generalized 
formulation.  Only  the  size  of  computer  core  and  time  requirements 
limit  the  detail  to  which  v  and  its  spatial  derivatives  can  be  represented 
in  the  model.  To  study  subduction  zones,  such  as  the  Kurile -Kamchatka 
arc,  it  is  convenient  to  represent  the  structure  with  a  two-dimensional 
grid,  norrml  to  the  arc,  in  which  velocities  are  specified  at  discrete 
points.  In  this  system  any  degree  of  structural  complexity  may  be 
represented,  limited  only  by  the  fineness  of  intervals  between  grid  points. 
The  wave  speed  and  its  derivatives  may  then  oe  determined  numerically 
at  any  point  within  the  structure  by  using  cubic  splines. 

It  has  been  instructive  to  study  a  short  slab  like  the  Aleutian  arc. 
which  has  reasonably  well-known  geometric  parameters,  before  proce¬ 
eding  to  the  longer  slab  of  the  Kurile  -  Kamchatka  arc.  Moreover,  the 
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existence  of  good  data  from  a  local  network  in  the  Central  Aleutians 
provides  a  unique  opportunity  to  study  a  range  of  possible  models  and 
plate  effects  that  have  application  elsewhere. 

Fig.  1.  1  is  a  cross-section  of  thirty-nine  intermediate  depth 
earthquakes  of  magnitude  4  or  greater  which  have  been  located  usim* 

O 

only  data  from  a  local  network  centered  on  Amchitka  Island  in  the 
Central  Aleutians  and  with  a  Herrin  spherically  symmetric  Earth  model. 
A  hypothetical  plate  model,  shown  for  reference,  was  independently 
derived  from  observed  travel  times  and  amplitudes  from  nuclear 
explosions  on  Amchitka  Island2’  3.  The  scattering  of  the  v/pocenters 
serves  to  illustrate  uncorrected  plate  effects  in  the  local  data.  The 
solid  circles  are  two  events  within  the  set  -  hich  weie  also  relocated 
but  with  a  common  network  of  twenty-nine  local  and  teleseismic  stations. 
In  addition,  an  attempt  was  made  to  correct  for  plate  effects  by  adjusting 
observed  arrival  times  by  their  corresponding  residuals  from  the 
Amchitka  explosions.  A  significant  reduction  in  the  standard  error  of 
an  observation  was  achieved. 


An  obvious  next  step  is  to  incorporate  the  plate  model  in  the 
hypo  center  solution  using  seismic  ray  tracing  in  the  manner  described 
by  En&dahl  .  In  Fig.  1.  2  the  thirty-nine  *vents  have  been  relocated 
using  a  relatively  simple  plate  representation  which  incorporates  a  7 
per  cent  increase  in  velocity  within  the  plate  over  Hie  surrounding  Herrin 
mn.'lc.  Th-  approximate  grid  size  for  this  model  was  5X5  km.  The 


Figure 
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effect  on  the  hypocenters  is  dramatic.  They  are  now  concentrated  in 
a  rather  thin  zone  (~  10  km)  in  the  upper  half  of  the  plate  which  also 
corresponds  to  the  colder,  more  brittle  region  of  the  slab  suggested 
by  thermal  anomalies  . 

A  more  realistic  velocity  model  may  also  be  constructed  from 
numerically  calculated  thermal  models  of  slabs.  In  Fig.  1.  3  is  shown 

3  Pl3le  mod°l  (on  a  10  X  10  »<nt  grid)  derived  by  perturbing  the  Herrin 
velocities  with  a  thermal  model  of  the  Aleutian  plate  derived  by  Sleep3. 
Hypocenters  relocated  using  the  local  network  and  this  model  again 
define  a  thin  zone  in  the  colder  upper  half  of  the  plate,  but  also  cor¬ 
responding  to  the  region  of  largest  velocity  (or  thermal)  gradient.  The 
solid  circles  are  the  same  events  shown  in  Fig.  1.  1,  except  relocated  by 
seismic  ray  tracing  using  this  model.  In  this  case,  station  corrections 
applied  Lo  the  observed  arrival  times  were  the  residuals  obtained  from 
the  Amchitka  explosions  by  tracing  through  this  model.  The  locations 
of  these  two  events  shown  in  Figs.  1.  1  and  1.3  are  in  close  agreement. 

This  is  to  be  expected  if  the  plate  structure  is  a  good  first  approximation 
to  the  real  Earth. 

The  tele-seismic  events  shown  in  Figs.  1.  1  and  1.  3  atso  provide  a 
good  opportunity  to  study  plate  effects  on  focal  mechanisms  and  dr/dA 
measurements  beyond  30’.  Figs.  l..|  and  1.5  are  comparisons  of  ray 
projeCons  on  the  lower  half  of  the  focal  sphere  (equal-area  azimuthal) 
for  the  spherically  symmetric  and  laterally  varying  Earth  models.  Only 
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igure 
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data  road  from  «  WKSS  stations  and  the  Amchitka  network  by  the 
investigate-  are  shown.  Changes  in  i  ,„d  =t  of  up  to  11’ and  20=  ,  respec¬ 
tively,  between  the  two  models  were  computed  for  rays  beyond  30\  The 
largest  changes  occurred  for  rays  departing  nearly  along  ,he  strike  of 
the  arc.  Th.s  can  be  seen  in  the  general  counterclockwise  rotation  of 
ray  points  to  North  America  in  the  NE  quadrant.  In  some  instances  a 
rotatton  of  focal  angles  of  this  magnitude  could  produce  a  corresponding 
change  in  the  strike  of  a  nearby  nodal  plane.  For  these  two  events  the 

distribution  of  ray  points  is  such  that  little  or  no  change  in  the  focal 
mechanisms  is  demanded. 

The  theoretical  dT/dAand  back  ar.imuth  of  the  wavefront  were 
also  computed  for  stations  beyond  30=  .  The  larges,  departures  from  a 
spherically  symmetric  model  were  0.02  sec/=  and  0.  4=  ,  respectively. 
This  is  no,  a  wholly  unexpected  result  since  only  a  small  par,  of  the 
telescismic  path  actually  secs  the  laterally  varying  mantle  at  the  source. 
Tits  does,  however,  suggest  quantitatively  that  plate  effect  will  not 
normally  be  detectabie  by  seismic  arrays  at  teleseismic  distances. 

g fatefully  acknowledge  the  cooperation  of  „,c  National  Center 
for  Atmospheric  Research  permitting  the  use  of  thoir  CDC  7600 
compute r  for  the  purpose  of  expediting  the  develop, non,  of  the  eompu- 
tational  techniques  used  in  this  research. 
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— — Staton  Corrections  for  Soviet  Stations 
E.R.  Tngdahl  and  R.  .^grawal 


Approximately  1000  residuals  of  lel^oismu-  P.nvos  reported  by 


•stations  in  the  Soviet  bulletins  for  earthquakes  deeper  then  300  km 


during  I  u70-  1(J73  were  examined  for  source -station  anomalies.  Resid¬ 


uals  were  first  separated  into  source  regions  by  station  to  determine 


if  mean  residuals  vary  with  the  body-wave  magnitude  of  the  earthquake 


(or  equivalently,  the  number  of  stations  reporting)  or  the  depth  for  each 


source  region.  No  obvious  effects  were  observed  in  the  plotted  data. 


at  least  in  part  because  the  Soviet  bulletins  are  selective  only  of  the 


larger  magnitude  earthquakes, 


Station  quality  appears  to  vary  roughly  with  reported  station 


sensitivities.  Stations  reporting  to  0.  1  sec  were  not  always  as  reliable 


as  those  reporting  to  0.  5  sec.  One  station  (MAG)  shows  an  abrupt 


change  of  mean  residual  of  about  3  sac.  during  April,  1971,  an  effect 


that  could  be  due  to  an  unreported  change  of  location. 


Source  biases  introduced  by  different  source  regions  were  removed 


by  computing  station  residuals  relative  to  a  reliable  base  station  (TLG). 


These  relative  residuals  had  means  consistent  within  a  standard  error  of 


0.  I  sec.  The  mean  residual  was  then  computed  over  all  source  regions 


for  the  base  station  and  the  relative  station  means  (to  the  base  station) 


adjusted  accordingly  to  obtain  absolute  station  corrections  relative  to 


the  Jeffrey  s-Bullen  tables.  The  results  are  plotted  in  Fig.  \ . 


The  bracketed  numbers  in  Fig.  2.  i  are  mean  station  corrections 


independently  determined  by  Kogan’  from  Eniwet  A  explosions  and 
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adjusted  by  +2.  0  sec.  to  the  J-B  tables.  .Although  these  are  small  samples, 
they  represent  high-quality  data  and  are  in  fair  agreement  with  the  results 
of  the  deep-earthquake  study.  In  the  region  north  of  the  Hindu  Kush  the 
agreement  is  to  within  0.  2  sec. 

This  study  demonstrates  the  usefulness  of  the  technique,  and  plans 
are  being  made  to  e.v*ond  the  results  to  other  Asian  stations  using  ISC 
data  from  1964-1970.  The  station  corrections  will  then  provide  a  means 
of  separating  out  purely  source  effects  from  observations  of  earthquaKes 
in  the  Hindu  Kush  and  Kurile -Kamchatka  regions. 

Rch  rence. 

1*  Kogan,  S.D.  (1960).  Travel  times  of  longitudinal  and  transverse 

waves,  calculated  from  data  on  nuclear  explosions  made  in  the  region 
of  the  Mr  rshall  Islands,  Izv.  ,  Geophys.  Ser.  ,  371-380. 
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3.  Depth  Determination  from  Long  Period  Body  Wave  Spectra 

G.  3oucher 

The  determination  of  the  depth  of  focus  of  earthquakes  is  one  of 
the  classical  problems  in  se. sinology.  Although  the  problem  may  be 
regarded  as  solved  in  principle,  through  use  of  body-wave  travel  times, 
depth  phases,  surface  wave  amplitudes  and  spectra,  and  so  forth,  well- 
known  practical  difficulties  remain,  particularly  in  the  case  of  shallow 
earthquakes.  In  this  research  Fourier  spectra  of  body  waves  recorded 
or.  standard  WWSSN  loi  -  period  seismograms  are  interpreted  in  terms 
of  the  spectral  modulations  produced  by  the  interference  of  near-source 
surface  reflections  and  the  direct  body  wave,  such  as  pP  and  P.  The 
method  of  analysis  has  practical  value  not  only  in  determination  of  depth 
of  focus  (which  may  be  useful  in  deciding  whether  a  suspicious  event 
might  be  an  explosion  or  an  earthquake),  but  also  in  determination  of 
the  focal  mechanism.  This  work  was  intended  to  lay  the  groundw'ork  for 
applying  the  process  to  digital  seismograms  from  the  High  Gain  Long 
Period  Network  (HGLt-'l.  It  has  been  found  that  a  relatively  simple 
model  of  a  body  wave  phase  and  one  or  more  near-source  surface 
reflections  produces  credible  and  consistent  results  in  some  cases, 
when  the  precise  timing  of  the  depth  phase  would  normally  be  in  doubt. 
Considerable  attention  has  been  given  to  maximization  of  the  chances 
for  success  of  the  method  and  to  the  sources  of  ambiguity  and  difficulty. 

Focal  depth  determination  by  means  of  long  period  body  waves 
complements  other  methods  of  focal  dep  h  determination.  When  only 
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teleseismic  data  are  available  to  define  the  soiree  parameters  of  an 
earthquake,  unless  the  depth  of  focus  is  great  enough  to  produce  unam¬ 
biguous  depth  phases,  determination  of  focal  depth  is  limited  to  indirect 
means.  Ba sically  this  is  because  of  the  trade -off  betu  een  focal  depth 

and  origin  time  in  teleseismic  location  algorithms.  Tsai  and  Aki^ap- 

0 

4* 

plied  a  theory  of  Sa ito  to  determination  of  focal  parameters  using 
surface  wave  spectra,  but  the  results  are  not  without  ambiguity,  and 
either  multipath  effects  or  multiple  events  at  the  same  source  may  give 
misleading  results.  The  problem  in  separating  depth-diagnostic  phases 
such  as  pP  from  P  for  shallow  earthquakes  is  that  the  seismic  signals 
are  too  close  together  in  time  to  be  distinguished.  With  short  period 
seismograms  the  problem  is  that  the  complicated  structure  of  the  earth 
results  in  a  body  wave  signal  that  lasts  for  several  tens  of  seconds. 
Because  of  the  narrow  bandwidth  of  typical  instruments  it  is  difficult 
to  pick  later  arrivals  from  the  nearly  sinusoidal  signal.  On  long  period 
records,  such  as  the  V  WSSN  long  period  instruments,  the  time  resolution 
is  relatively  poor,  and  the  signals  are  somewhat  smeared  out  in  time 
because  of  the  amplitude  and  phase  response  characteristics  of  the 
instruments,  although  the  s; goals  do  not  appear  so  complex.  Thus  even 
if  it  is  apparent  that  a  mixture  of  phases  is  present  in  a  packet  of  body 
waves,  it  is  quite  another  matter  to  identify  their  arrival  times  and 
relative  senses  with  certainly.  It  has  been  learned  in  this  study  that 
the  application  of  Fourier  spectrum  techniques  to  long  period  body- 
waves  of  earthquakes  shallower  than  about  SO  km  makes  the  focal  depth 
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problem  tractable.  Part  of  the  following  discussion  will  be  devoted  to 
exploring  the  liberties  that  may  be  taken  with  mathematical  rigor  and 
physical  reality  in  order  to  achieve  useful  results. 

After  treatment  of  the  simple  theory,  a  few  examples  will  be 
shown  to  demonstrate  the  expected  modulations  on  the  body-wave 
spectrum.  A  second  approach  will  be  a  comparison  of  real  seismograms 
with  synthetic  seismograms  made  by  convolving  theoretical  seismic 
pulses  and  their  echoes  with  the  instrument  response.  The  focal  depth 
and  corner  frequency  may  be  determined  simultaneously  by  visual  trial 
and  error,  without  recourse  to  Fourier  analysis. 

Theory. 


An  easy  way  to  conceptualize  the  use  of  spectrum  modulations  to 
characterize  secondary  arrivals  is  to  consider  the  convolution  of  a  signal, 
say  an  elementary  P  wave,  with  a  pair  of  impulses  which  may  or  may  not 
be  mutually  reversed  in  sense.  It  will  be  assumed  that  for  practical 
p,  noses  the  signals  to  be  separated  are  similar  in  shape,  but  with  pos¬ 
sibly  different  amplitudes.  The  amplitude  spectrum  of  a  pair  of  identical 
impulses  is  a  sinusoidal  curve  with  either  a  maximum  or  a  minimum  at 
zero  frequency,  depending  upon  whether  the  impulses  are  of  the  same  or 
opposite  sense,  respectively.  In  the  general  case  where  the  relative 
p  ‘  'plitude  of  the  second  signal  is  +  a,  and  the  Fourier  amplitude  spectrum 
of  the  common  pulse  is  f  (a),  where  *  is  angular  frequency,  then  the 
modulated  Fourier  spectrum  due  to  the  interference  of  the  second  arrival 
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is  f  (i)  .  (]  ta  +  2a  cos  (XT)  where  T  is  the  time  interval  between 
events  ’  .  If  the  two  arrivals  are  dissimilar  the  situation  in  the 

frequency  domain  becomes  much  more  complicated.  For  moderate 
earthquakes  such  as  those  studied  here  the  soi  e.  dimensions  are  small 
enough  that  the  separate  pulse  arrivals  may  be  regarded  as  similar  for 
practical  purposes. 

The  utility  of  the  method  is  aided  by  the  form  o'  ‘  elementary 
source  spectrum  of  the  earthquake.  Common  to  dislocation  models  of 
the  earthquake  source  function  is  an  essentially  flat  spectrum  from  some 
'corner  frequency'  down  to  zero  frequency  ’  .  Thus  modulations  due  to 

interference  predominate  over  expected  peculiarities  of  the  source 
mechanism  as  long  as  the  frequency  range  of  interest  is  below  the  corner 
frequency.  Because  of  the  shape  of  the  instrument  response,  noise,  and 
a  technical  difficulty  pointed  out  by  Lin  He  and  Sacks',  it  is  difficult  to 
establish  the  very  low  frequency  points  of  the  spectrum.  On  the  other 
hand,  at  sufficiently  high  frequencies  attenuation  and  the  effects  of  the 
crustal  response  become  severe.  Thus  there  is  a  region  of  the  spectrum 
where  the  frequencies  are  neither  too  high  or  too  low,  where  path  and 
multiple  phase  effects  predominate.  One  goal  of  this  study  was  to  de¬ 
termine  the  region  of  the  spectrum  where  this  is  true  for  shallow  events 
of  moderate  size  as  recorded  by  \\  WSSN  long  period  seismographs. 

In  the  case  of  long-period  seismograms,  a  recorded  body-wave, 
often  can  be  represented  b)  the  convolution  of  the  instrument  response 
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with  a  rather  simple  pulse  shape  ot  shirt  duration.  1  hat  is.  the  effect 
ol  layering  or  other  i nhonmgeueil  ies  alon;4  (lie*  p ropaitut ion  path  is 
sometimes  negligible,  so  that  the  effect  of  removing  the  instrument 
response  is  to  reproduce  essentially  the  form  of  the  seismic  pulse. 

Thus  it  may  be  effective  to  remove  the  instrument,  response  from  the 
complex  Fourier  spectrum  of  the  signal.  The  reason  for  this  ii  that 
the  smearing-out  of  the  seismic  signal  by  the  instrument  response  is 
then  eliminated,  and  seismic  pulses  of  short  di  ration  c  ictually  be 
'read'  from  the  deconvolved  seismogram  of  direct  ground  motion.  This 
is  usually  less  satisfactory  than  looking  at  the  spectrum  to  determine 
delay  and  sense  of  multiples  of  a  pulse,  but  it  is  instructive  in  other 
ways  since  the  phase  information  is  utilized,  giving  a  more  complete 
picture.  It  can  also  be  diagnostic  of  cases  where  the  simple  theory  is 
inadequate,  pa rticuln  rly  when  extra  pulses  are  present. 

although  problems  of  noise  and  multiple  events  are  discussed 
in  the  sections  on  experimental  results,  a  few  comments  on  the  practical 
aspects  of  using  Fourier  transforms  for  non-stationary  signals  in  the 
presence  of  noise  may  be  in  order.  The  essential  point  is  that  a  seis¬ 
mic  body  wave  signal  or  phase  is  finite  in  duration.  On  a  raw  digitized 
seismogram,  if  the  duration  of  a  source  pulse  is  only  a  few  seconds, 
then  a  data  sample  no  longer  than  20  or  30  seconds,  which  is  about  the 
length  of  the  most  important  part  of  the  impulse  response  of  the  seis- 
mograoh  system,  may  be  adequate.  To  take  a  longer  sample  is  simply 


-  1 1 


to  include  other  possible  phases,  and  more  significantly  to  include  a 
greater  proportion  of  noise  energy.  It  is  obvious,  of  course,  that  the 
sample  must  be  long  enough  to  include  a  sufficient  sample  of  all  the  pulses 
of  interest.  For  example,  a  sample  that  includes  an  adequate  length  of 
P,  but  only  the  beginning  of  pP  will  be  modulated  correctly  in  the  upper 
frequency  range,  but  he  bottom  end  of  the  spectrum  will  be  incorrect, 
and  the  vital  first  hole  may  be  missing.  But  inclusion  of  unwanted  phases 
produces  unwanted  modulation  of  the  spectrum,  making  interpretation 
difficult.  Since  microseismic  noise  is  of  a  quasi-stationary  character 
the  noise  power  in  the  spectrum  is  nearly  proportional  to  record  length. 

For  signal-to-noise  ratios  of  the  order  of  10,  as  is  the  case  for  earth¬ 
quakes  in  the  size  range  studied  here,  any  unnecessary  amount  of  noise 
is  intolerable.  For  isolated  signals  signal-to-noise  ratio  is  a  decreasing 
function  of  sample  length  beyond  some  point.  This  is  in  contrast  to  the 
situation  for  stationary  signals  and  noise  whore  the  signal  to  noise  ratio 
doesn't  vary  much  with  sample  length,  and  where  it  may  be  possible  to 
characterize  the  noise  well  by  developing  a  statistically  better  deter¬ 
mined  average.  In  the  case  of  seismograms,  the  noise  may  not  be  stationary, 
and  no  certain  advantage  accrues  from  studying  noise  samples  that  are 
longer  than  adequate  signal  samples.  This  is  especially  true  when  the 
record  contains  noise  trans  ents  or  interfering  phases,  as  often  turns 
out  to  be  the  case.  It  is  critical  to  bear  in  mind  that  the  travel-time 
curve  of  any  particular  phase  may  be  complicated  in  certain  distance 
ra  ges,  with  duplication  o r  triplication,  cusps  and  other  phenomena. 


some  of  which  may  not  be  understood.  For  this  reason  the  interpre¬ 
tation  of  an  earthquake  involves  numerous  lines  of  evidence,  im  .hiding 

as  many  phases  as  possible  and  as  many  stations  at  different  distances 
as  may  be  necessary. 

\\  ilh  reference  to  the  problem  of  discrimination  of  underground 
nuclear  explosions,  it  should  he  noted  that  in  the  case  of  P  waves  and 
other  waves  of  the  P  type,  there  will  be  very  little  energy  at  long  period 
because  of  the  cancellation  effect  of  the  ncar-source  surface  reflection 
of  the  P  type.  Thus,  an  interesting  special  case,  to  be  discussed  in 
Section  4,  is  that,  in  which  no  detectable  long  period  signal  is  present. 
Experimental  Data  and  Observations. 

Data  for  the  numerical  experiments  were  obtained  by  digitization 
of  photographic  enlargements  made  from  WWSSN  film  chips,  using  a 
Calma  digitizer.  The  digitizer  is  claimed  to  have  accuracy  of  0.003 
inch;  the  principal  limitation  upon  the  accuracy  of  the  digitization  is 
the  steadiness  of  the  operator.  The  digitizing  rate  was  always  at  least 
1  sample  per  second  of  record  time,  and  usually  4  per  second.  Aliasing 
effects  appear  to  be  negligible.  Below  0.  5  Hz  the  spectra  for  records 
digitized  at  4  samples  per  second  were  not  significantly  difficult  in  the 
spectrum  below  0.  5  Hz  from  those  with  a  sampling  rate  of  1  per  second. 
The  Fast  Fourier  Transform  algorithm  was  used  to  compute  Fourier 
spectra,  using  either  128or  512  points,  and  filling  in  with  zeroes  if  the 
sample  was  shorter  than  the  specified  number  of  transformed  points. 


This  study  was  based  on  earthquakes  located  in  the  central 
Aleutian  Islands,  principally  because  data  from  the  Aleutian  Island 
Seismic  Network  helped  to  insure  that  the  depths  obtained  from  tele- 
seismic  solutions  are  fairly  accurate.  In  this  way  the  technique  of 
depth  determination  from  long  period  body  waves  can  be  evaluated. 

Data  from  stations  as  near  as  21  degrees  (COL)  and  as  far  as  90 
degrees  (SJG)  were  used.  Magnitudes  of  events  were  in  the  range  5.  5 
to  6  as  determined  by  the  National  Earthquake  Information  Service  of 
the  U.  S.  Geological  Survey.  The  depths  given  in  the  Preliminary 
Determination  of  Epicenters  were  used  as  a  guide  in  preliminary 
selection  of  data.  It  is  mainly  by  comparison  with  the  depths  deter¬ 
mined  by  the  PDE  program  that  the  method  is  justified  for  use  in  other 
areas  where  depth  control  is  poorer. 

TABLE  I 


Event 

Lat. 

Long. 

Depth, km 

m/M 
b  s 

Location  Figu 

05-14-69 

19  32  54. 2 

51.  3N 

179.  2AV 

21 

6.  2/7.  0 

Andreanof  Is. 

09-12-69 

08  57  07. 3 

51. 2N 

179. 2  W 

48 

6.  0/6.  6 

1  1 

07  - 18-70 

01  48  38.9 

51. 4N 

1 7  S .  5\V 

46 
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1  1 

08-04-69 

17  19  19.6 

5.  7S 

125.  3E 

521 

m.  -6.2 
b 

Banda  Sea 

08  -  30 -70 

i  *7  i  /I  n 
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151.6E 

64  5 

m. -6. 6 
b 

Sea  of  Okhotsk 

OH  -  30 -70 
17  4o  09. 0 
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Parameter,  for  the  earthquake,  treated  i„  this  report  are  given 
a  Jlc  I.  Tne  main  emphasis  is  upon  two  groups  of  earthquakes, 
rather  shallow  (20-25  km)  and  somew hat  dceper(40-60  km).  Brief 
mention  will  be  made  of  a  deep  focus  event,  net  to  study  pP,  but  to 

illustra te  some  apparent  source  mechanism  effects  that  are  susceptible 
to  the  same  sort  of  analysis. 

The  theory  of  detecting  PP-P  or  sP-P  intervals  from  modulations 
of  the  spectrum  is  straightforward .  but  typical  spectrums  made  from 
seismog rams  have  mediations  from  other  causes  that  must  he  sorted 
out  before  the  correct  ones  con  be  identified.  The  spectrum  will  have 
irregularities  due  to  background  noise,  effects  o,  the  transmission  path 
«ha.  were  neglected  in  the  theory,  and  perhaps  unsuspected  multiplicity 
oi  the  source.  The  s P  phase  is  also  of  primary  importance.  The  focal 
mechanism  of  .he  earthquake,  which  affects  both  the  relative  amp.ffude 
and  sense  of  the  phases  being  examined,  adds  another  set  of  variables 
.ha.  affect  each  setsmogram  differently.  The  solution  to  the  problem 
g  the  cornea  set  of  modulation,  is  t„  „sc  a  sot  of  different 
sample  lengths  and  time  win*,,-,  for  each  recorded  phase  (chosen  logically 

'r0“‘  lh°  0r‘*“U"  *0t,m0*ram-  »'  «"•••-•>.  •'"‘1  to  compare  observations 

am°nt!  *  °'  r"0‘'d“'^  f'°"'  stations.  If  available,  a 

nearby  event  can  serve  as  a  control.  The  ftrs,  technique  eliminates 

spurious  modulations  that  arise  somehow  from  .he  choice  of  sample 

Jcneth.  for  example  when  an  undesired  phase  is  included  within  ,he  sample. 

The  second,  to  use  a  specific  example,  exploits  the  small  but  measurable 
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variation  of  pP-P  interval  with  epicent ml  distance  as  predicted  by 
the  ti  ivcl  time  tables,  t.learly  a  hypothetical  choice  of  a  particular 
spectral  modulation  that  is  not  only  stable  with  respect  to  sampling 
length,  but  also  shows  the  proper  behavior  as  a  function  of  epicentral 
distance  is  likely  to  be  the  correct  one.  A  satisfying  exercise,  in.  fact, 
is  to  vary  the  length  of  the  sampling  window  about  an  appropriate  phase 
^  order  to  deliberately  include  an  unwanted  phase.  As  the  window’ 
length  is  extended  to  include  the  intruding  phase,  the  effect  of  that 
phase  can  be  seen  to  develop  in  the  spectra.  This  serves  as  encour¬ 
agement  that  the  technique  is  basically  correct.  It  also  shows  that,  as 
long  as  one  does  not  stray  from  the  frequency  range  where  the  data  are 
reliable,  the  sampling  problem  is  not  too  severe, and  it  is  possible  to 
use  rather  subjective  criteria  to  determine  just  what  constitutes  the  wave 
packet  of  interest,  and  what  should  be  rejected.  One  reason  this  point 
must  be  emphasized  is  that  for  the  earthquakes  of  primary  interest, 
the  signal  to  noise  ratios  are  typically  rather  low.  The  second  is  that, 
for  earthquakes  in  the  upper  end  of  the  magnitude  range  studied,  i.  e. 
greater  than  6.  0,  the  likelihood  that  the  event  will  be  a  muUiple  one  is 
considerable.  1  hat  situation  must  be  identified  by  inspection  of  several 
seismograms  or  else  from  indications  given  by  the  spectra.  It  is  neces¬ 
sary  to  look  at  records  from  several  stations  and  from  several  different 
phases,  if  possible,  to  resolve  the  usual  ambiguities  that  are  seen. 

The  first  example  considered  is  the  rather  complicated  earth- 
qua.  e  of  May  14,  1969.  The  whole  suite  of  spectrum  data  obtained  are 
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.horn,,  in  order  to  provide  a  complete  picture  of  the  interpretive  process. 
Other  examples  trill  then  be  discussed  in  less  detail,  representing  pre - 
suma bl y  less  complex  situations:  Figure  3.  1  shorts  the  entire  suite  of 
data  from  stations  COL  (P).  KIP  (P,  PP;  COL  (P).  and  SJG  (P).  The 
letters  in  parentheses  here  signify  the  phase  which  together  with  its 
surface  reflection,  was  used  in  the  interpretation:  e.  g.  ,P)  means  the  pP-P 
(or  possibly  sP-P)  time  interval  was  determined.  Figures  3.  2  and  3.  3 
illustrate  other  examples  as  described  in  the  figure  captions. 

It  is  sometimes  unavoidable  to  include  unwanted  phase  arrivals 
within  the  time  window  used  to  determine  the  spectrum.  Since  the  sP 
phase  is  likely  to  be  present  in  addition  to  pP,  for  example,  the  spectrum 
is  sometimes  modulated  by  more  than  1  sinusoid.  Interpretation  by 
inspection  of  a  sum  of  sinusoids  of  mixed  phase  and  different  amplitudes 
is  difficult,  so  that  a  useful  approach  is  to  model  various  interpretations 
by  a  sequence  of  impulses  in  the  time  domain,  which  are  readily  converted 
to  spectra  for  comparison.  Since  the  instrument  response  and  the  source 
spectrum  in  the  period  range  of  10  to  100  seconds  are  reasonably  flat, 
it  does  not  introduce  much  error  to  approximate  phase  arrivals  by  spikes 
for  moderately  large  earthquakes,  the  main  proviso  being  .ha.  the  corner 
frequency  is  above  the  frequency  range  of  interest. 

inasmuch  as  a  some  subjective  analysis  turns  out  to  l.e  required 
tn  most  cases  in  the  interpretation  of  spectra,  and  since  the  original 
seismograms  are  the  primary  resource  on  which  to  base  the  subjective 
analysis,  an  alternative  method  of  determining  depth  of  focus  suggests 
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Fig.  3.1.  Earthquake  of  'lav  14,  1969.  Spectra  of  teleseismic  P  waves  with  dif¬ 
ferent  sample  lengths  to  show  the  existence  of  a  stable  minimum  (dot)  at  0.133- 
0.1*1  Hz,  signifying  a  pp  -  P  • ime  between  7.1  and  7.5  sec.  COL  =  College .Alaska , 

KIP  -  KL|  Hawaii;  GOl.  =  Golden,  Colorado;  SJC  =  San  Juan,  Puerto  Rico.  The 

choice  of  pP-F  time  is  also  based  upon  additional  phases  not  show*:  here,  sim  i  as 
pPP-Pl’  and  ^J-S.  Spectra  are  normalized  to  maximum  apmlituue,  shown  by  am  s. 
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ig.  3.2.  Earthquake  of  Sept.  1?.,  1909.  Spectra  of  P  waves  showing  pP  out  of  phase 
vi.'h  p  at  Matushiro  (MAT)  and  College  (COL),  and  in  phase  at  Kipapa  (KIP),  The 
mi'.ihers  give  sample  lengths  in  seconds.  The  short  sample  at  KIP  was  inadequate  to 
develop  the  low  frequency  hole  corresponding  to  the  first  minimum.  The  pP-P  time 
estimated  from  these  and  other  data  is  9.5-10  sec. 
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Hg.  3.3.  Earthquake  of  July  18,  1970.  Teleseismic  P-wave  spectra.  A  and  B  are 

!r,"Vn'  Golden,  Colorado  record.  Sample  lengths  are  13  and  67  sec,  respectively 
C  ,mJ  D  are  from  the  College,  Alaska,  record.  Sample  lengths  are  18  and  43  sec 
respectively .  This  earthquake  was  located  by  the  National  Earthquake  Information 
["'  “I  f  k«n  depth.  The  sample  length  in  the  short  samples  was  inadequate 
to  .  i splay  ..ay  but  1  of  the  m.-.m  associated  with  the  probable  pi’-p  -ime 
which  is  rather  well  determined  at  11.3  sec  from  the  long  sample  at  Colle-e 
and  17.7  sec  at  Golden.  °  * 
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itSC“'  T'US  "  PrCPare  «‘-mogr.ra.  based  upon  different 

as  sumptions  about  the  relative  sire,  timing  and  sense  o,  arrival  o, 

Phases,  using  a  reasonable  model  for  the  pulse  shapes  to  be  expected. 

The  sequence  of  pulses  is  then  convolved  with  the  instrument  response 

and  plotted  to  produce  an  overlay  for  direct  visual  comparison  with  the 

seismog  ram.  ,n  practice,  i,  turns  on,  to  be  convenient  and  necessary 

inCl"dC'  ^  "*  PH..e  as  par,  of  the  process,  sine 

WWSSN  l0  'fi  PeH°d  il’StrUme!’t“  *»  frequencies  of  the  order 

°f  1  Hr.,  even  though  a,  reduced  magnification.  Thus  a  rather  simple 

process  can  be  used  to  determine  not  only  depth  of  focus  but  corner 

frequency  in  favorable  cases.  Figure  i  a 

°  C  3,4  shows  3  representative  set  of 

synthetic  seismograms  for  Hiffar-nni 

■  tor  different  assumptions  using  the  far -field 

pulse  shape  suggested  by  Brun_e‘.  For  „se  in  other  phases  ^  ^ 

”  mier°mm  Pl°‘  °f  -H  artifical  seismograms  was 

prepared  at  a  very  sm.,11  cost,  clearly  it  is  not  possible  to  cover  all 
possible  situations  with  a  finite  number  of  overlay  seismograms,  but 
a  Pood  estimate  of  earthquake  parameters  can  be  made  with  a  fairly 
coarse  spacing  of  sample  parameters.  The  time  of  first  aero  crossing, 
th=  ratio  of  the  overshoot  pulse  amplitude  (second  half  cycle)  to  the  first 
half  cycle  an, p, Unde,  and  the  ,ime  of  the  second  r.cro  crossing  are  the 
haste  Parameters  used  to  make  an  interpretation.  Furthermore, 

.noma  1  ,us  arrivals  are  some, hues  . . .  „  a  glance  from  the  com- 

pansnn.  To  ill,, s, rate  the  usefulness  of  the  technique.  Fig.  5.  5  shows 
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3.5.  Selected  P-wavos  from  the  earthquake  of  Aug.  4,  1969,  showing  the  possibility  of  f:tt'>-  -n'Setic  puls 
ith  different  corner  frequencies.  Here  c.  replaces  the  K  in  figure  4,  and  is  roughly  equal  to  the  reciprocal 
rner  frequency.  Since  this  is  a  deep  earthquake  (645  km)  no  surface  ref lecticr.s  like  p?  or  sP  are  incl".  -d  « 
-pie  shown,  yost  theories  predict  that  the  corner  frequency  will  be  slightly  different  _n  difrerent  diructi' 
on  the  source. 
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an  example  of  a  deep  earthquake  in  the  Banda  Sea  with  simple  pulse 
shapes  In  a  second  earthquake,  Fig.  3.6,  a  secondary  phase 
apparently  opposite  in  sign  to  the  initial  P  phase  is  observed  at  about 
1  1  seconds  delay.  The  phase  is  observed  on  several  stations  with  con¬ 
sistent  timing,  so  that  it  does  not  appear  to  be  a  result  of  some  local 
peculiarity  of  t lie  travel  path.  The  interpretation  of  this  phase  will  not 
be  considered  here. 

Conclus  ions . 

This  first  conclusion  of  this  study  is  the  general  success  of  the 
method  as  applied  to  standard  \\\ .  SSN  long  period  seismograms  and  the 
determination  of  its  range  of  applicability.  Second  is  the  identification 
of  the  sources  of  probable  error,  and  methods  of  dealing  with  them.  Third 
is  the  recognition  of  directions  in  which  to  seek  improvement  of  the  limiting 
parameters  to  make  the  method  more  generally  useful. 

The  principal  conclusion  from  the  study  is  that,  given  certain 
preconditions  on  the  events  themselves,  it  is  possible  to  identify  specific 
modulations  in  earthquake  body-wave  spectra  obtained  from  digitized 
paper  records  of  WVSSN  long  period  instruments.  In  general  it  is  safe 
to  say  that  this  can  hr  done  in  the  range  of  frequencies  where  the  signal- 
to -noise  ratio  is  ot  the  order  of  10  or  more,  where  noise  must  be  con¬ 
sidered  to  arise  both  from  (lie  earth-instrument  system  on  the  one  hand, 
and  from  the  process  ol  digilirat !<"  on  the  other.  Taking  a  conservative 
but  realistic  position  that  the  useful  frequency  range  is  0.01  to  0.2  Hz, 
depth  as  shallow  as  6  or  10  km  from  sS-S  and  10  or  15  km  from  pP-P, 


depending  respectively  t.pon  whether  the  surface  reflection  is  in-phase 
rhssc  with  the  direct  wave. etui  be  determined.  Absence  of 
any  modulation  may  be  useful  information,  but  is  no.  definitive  by  itself. 
Sometimes  only  a  single  'hole'  can  be  located  within  the  reliable  frequency 
range.  The  interpretation  of  the  spectrum  may  be  aided  by  a  determi¬ 
nation  of  the  focal  mechanism  of  the  earthquake,  princ.pally  to  determine 
whether  r.ero  frequency  is  expected  to  be  a  maximum  c  r  a  minimum. 
Conversely,  the  spectral  modulation,  if  determined  veil  enough,  may 
contribute  to  determination  of  the  focal  mechanism,  by  indicating  the 
firs,  motion  of  various  phases  emerging  front  the  top  of  the  focal  sphere. 

A  technique  that  deals  with  noise  problems  about  as  well  as  careful  analysis 
of  the  spectra,  and  takes  advantage  of  some  of  the  useful  integrating 
properties  of  the  human  eye.  ts  compare  a  phase  on  the  original  seis¬ 
mogram  with  a  synthetic  pulse  generated  from  a  pulse  and  suitably  timed 
echo  of  the  proper  amplitude.  This  simple  procedure,  which  by  a  natural 
extension  includes  an  estimation  of  the  corner  frequency  of  the  event, 

produces  surprisingly  good  results  and  might  well  become  a  useful  par, 

Of  routine  interpretation  of  seismograms. 

The  discussion  of  the  sources  of  probable  error  will  make  it  ap¬ 
parent  why  this  particular  technique  has  rarely,  if  ever,  been  applied  in 
the  past.  Such  problems  as  micros eismic  noise,  instrument  problems 
(such  as  spurious  transients  I,  and  errors  in  digitizing  paper  records  are 

haxe  piedictabU*  results  on  the  spectrum.  More  fundamen¬ 
tal  problems  arise  from  properties  of  the  earthquake  source  and  the  earth 
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itself.  The  most  significant  manifestation  of  the  former  is  the  like!. hood 
that  an  "earthquake"  may  in  fact  consist  of  a  number  of  major  events 
which  can  be  sorted  out  only  by  considerable  effort.  In  the  present  study 
ol  shallow  to  intermediate  depth  ea  rthquakr-s  in  the  central  Aleutian 
Islands  region,  it  appears  that  events  with  magnitudes  (m^)  near  6.0  and 
above  ,irc  quite  likely  to  exhibit  multiplicity,  whereas  for  events  of  perhaps 
one  -half  magnitude  less,  multiplicity  seems  much  rarer.  Of  course,  if 
multiplicity  is  the  object  of  ttuly,  spectral  interpretation  techniques  can 
easily  be  applied,  using  appropriate  caution.  The  difficulty  presented  by 
the  Earth  itself  as  a  transmitting  medium  is  obvious  from  consideration 
of  a  single  seismogram.  From  the  viewpoint  of  classical  seismology, 
there  are  a  number  of  'arrivals'  visible  on  a  seismogram,  many  of  which 
can  be  identified  on  the  basis  of  classical  ray  theory,  or  as  leaking  modes 
or  surface  wave  modes.  Unfortunately  every  wiggle  on  every  seismogram 
has  not  been  accounted  for  in  detail,  so  that  there  may  be  uncertainty 
about  the  source  of  some  of  the  complications  visible  on  a  spectrum.  But 
even  avoiding  this  particular  complication,  the  presence  of  other,  well 
known  phases  on  the  record  near  the  one  of  interest  means  that  the  time 
sample  used  for  the  Fourier  transformation  must  be  of  limited  length,  and 
it  also  means  that  the  baseline  of  the  record  will  be  poorly  determined, 
lienee  the  spectral  amplitudes  at  very  low  frequencies  are  not  useful  as 
a  rule,  and  a  key  point  in  the  interpretation  is  always  the  presence  of  a 
minimum  or  a  maximum  .it  zero  frequency.  As  the  time  sample  is 
lengthened  to  include  more  and  more  ot  an  undcsiied  phase,  tie  effect 
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can  be  seen  on  the  spectrum  as  an  increased  degree  of  modulation  of 
its  own  which  eventually  may  swamp  the  modulation  that  is  being  sought 
for  interpretation.  The  practical  solution  to  all  of  these  problems  is  to 
compare  various  types  of  data.  For  each  record,  a  proper  sample  length 
or  range  of  sample  lengths  must  ^  “hosen.  i.i  order  to  eliminate  spurious 
arrivals  and  yet  have  as  long  a  sample  cf  data  as  is  practicable,  to  improve 
the  reliability  of  the  lower  frequency  ranges  of  the  spectrum.  It  is  further 
necessary  to  use  as  many  different  phases  as  possible,  and  as  many  records 
as  may  be  necessary  to  make  sure  the  interpretations  are  consistent  and 
reliable. 

Since  the  motivation  for  the  use  of  long-period  seismograms  in  the 
first  place  was  to  avoid  the  complications  of  the  short-period  crustal 
response  of  the  earth,  the  only  useful  extension  of  the  usable  frequency 
range  to  higher  frequencies  would  be  in  the  direction  of  suppressing  the 
mic roseismic  noise  of  periods  near  and  just  below  10  seconds,  since  that 
w’ould  allow  accurate  determination  of  the  spcctruins  of  events  of  smaller 
magnitude  and  lesser  depth  than  is  currently  the  case.  Some  improvement 
at  tile  longer  period  end  of  the  spectrum  would  presumably  result  from  use 
of  instruments  of  improved  stability  at  periods  greater  than  100  seconds, 
preferably  recorded  digitally  to  bypass  t lie  hand  digitization  process. 

Clearly  the  better  the  Jong  period  end  of  the  spectrum  can  be  specified, 
the  easier  it  is  to  determine  whether  zero  is  a  maximum  or  minimum,  and 
tiie  more  accurately  spectrum  modulations  can  be  identified.  A  case  of 
considerable  interest  is  that  of  underground  explosions.  Since  nearly 
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tm.nl  (Instructive  interference  a,  lone  P.  rind,  results  frem  the  interference 

energy  ,1  lone  periods  in  the  1' 

phase  of  an  event,  even  at  a  ample  station,  would  he  indicative  of  an 
explosion  or  at  least  an  extremely  shallow  event,  although  tectonic  strain 
release  could  of  course  negate  the  converse  proposition. 
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4'  iI^_oQlGI.P  Digit., 1  a,l.,  fo,-  „f  |ioilv  or  F.>rthq:uke< 


and  Explosions  in  Central  /sia 
G.  Boucher 


The  objective  of  this  program  is  to  explore  the  possibility  of  ap- 
plying  the  spectrum  modulation  techniques  discussed  in  Section  3,  or 
other  techniques  that  might  prove  fruitful,  to  digital  data  from  the  High 

Gain  Long  period  Seismograph  Network.  This  work  m  currently  in 
progress. 


In  principle,  long  period  body  waves  of  a  seismic  event  may  be 
expected  to  contain  useful  information  about  source  size,  focal  mechanism, 
and  depth  of  locus.  In  practice  there  is  a  lower  limit  to  the  magnitude  of 
an  event  lor  which  long  period  body  waves  can  be  observed,  and  in  general 
this  limit  is  larger  than  the  lower  limit  of  events  for  which  surface  waves 
can  be  detected,  at  lea st  for  shallow  earthquakes.  However,  because  of 
the  potential  usefulness  of  the  long  period  body  wave  information,  it  is 
fruitful  to  determine  the  circumstances  under  which  the  waves  can  be 
observed,  then  to  extract  the  information  that  is  available.  It  is  neces¬ 


sary  to  extend  the  study  to  the  smallest  possible  events,  in  the  context 
of  the  problem  of  discrimination  of  underground  nuclea-  explosions, 
because  discrimination  in  the  current  state  of  the  art  becomes  less  certain 


with  decreasing  magnitude  of  event.  There  are  two  theoretical  advantages 
inherent  in  tin  digital  data  from  the  I1GLP  network.  First,  the  instruments 


Have  been  designed  lor  increased  stability  in  the  intermediate  and  Ion- 
period  range.  Second,  the  digital  instruments  have  wide  dynamic  range 


and  the  data  are  convenient  for  v 


annus  computer  manipulations  of  the 


signal,  including  rotation  of  coordinates,  filtering,  and  Pnear  as  well 
ms  non-linear  operations  in  the  time  and  frequency  domains. 


In  the  initial  stage  of  this  investigation,  the  approach  may  be 
likened  to  that  of  the  classical  seismologist.  That  is,  records  from 
various  events  and  stations  are  examined  to  determine  the  presence  and 
the  characteristics  of  various  body  phases,  and  to  characterize  the  back¬ 
ground  noise.  One  of  the  striking  differences  between  the  presumed 
explosions  and  the  small  number  of  earthquakes  studied  so  far  is  their 
gi  cully  di  tic  rent  exeit.it  ion  of  long  period  energy  of  various  types,  using 
the  short -period  body  wave  mgnitude  m^  (as  determined  by  the  PDE 


calculations)  as  the  parameter.  This  is  m  accord  with  previous  studies, 

but  it  is  particularly  visually  striking  on  HGLP  seismograms.  For  example, 

an  earthquake  with  given  as  6.0  commonly  has  body  wave  amplitudes 

10  or  more  times  the  background  noise  amplitude  and  it  is  simple  to  nr  ve 

down  the  record  and  pick  out  phase  after  phase  ar  oredicted  by  the  travel 

time  curves.  A  typical  explosion,  on  the  other  hand,  with  m  given  as 

b 

6.0,  is  far  loss  impressive.  In  particular  phases  of  the  P  type,  including 
multiples  and  surface  reflections  generally  cannot  be  identified  at  all  on 
the  raw  records.  It  is  frequently  possible,  to  identify  a  few  phases  of 
the  SV  type,  notably  S  and  PS  (or  more  likely  SP)  in  appropriate  distance 
lunges,  but  they  arc  almost  lest  in  Ihe  noise  in  most  cases,  and  the 
multitude  of  phases  predicted  by  ill  travel-time  curves  are  generally 
nol  discernible,  /'ll  attempt  to  quantify  this  difference  will  be  made;  it 
is  considered  to  bo  a  significant  body  of  information  to  contribute  to  the 


problem  of  earthquake  -  explosion  d  is,  riminat  ion.  rhe  excitation  of 
long  period  waves  of  the  SY  type  by  explosions  is  considered  rather 
striking,  and  will  be  the  subject  of  a  special  commentary. 

The  research  program  to  date  has  been  devoted  prima rily  to 
analysis  of  presumed  explosions  in  Central  Asia.  Earthquakes  in  the 
general  region  have  also  been  studied  in  order  to  give  a  broader  frame¬ 
work.  This  initial  emphasis  has  necessarily  directed  the  bulk  of  the 
effort  into  certain  categories,  because  of  the  signal  amplitude  constraints 
imposed  by  the  existing  suite  of  explosion  data.  The  largest  recent 
explosions  within  the  Soviet  Union  have  been  assigned  magnitudes  (  ) 

of  6.  3  at  most.  Such  explosions  produce  rather  small  body  waves  on  the 
IIGLP  instruments  (with  respect  to  the  background  noise  levels)  and  so 
far  in  this  study  only  surface  waves  have  been  observed  from  explosions 
with  magnitudes  as  small  as  5.  0.  For  this  reason  the  thrust  of  the  analysis 
has  been  in  the  direction  of  identifying  and  comparing  different  wave 
arrivals  in  order  to  establish  general  familiarity  with  the  character  of 
the  seismograms;  spectrum  analysis  and  optimized  filtering  of  the  records 
will  be  the  subject  of  the  remainder  of  the  work  on  the  project. 


A  constraint  upon  the  progress  of  the  research  which  had  been 
unexpected  at  its  inception  has  been  the  relatively  small  amount  cf  data 


actually  available,  compared 
of  the  MG  UP  type  are  now  in  e 
the  author  has  ever  seen  for  a 
is  substantially  less  than  that. 


with  the  potential  amount.  Some  15  stations 
xislence ,  but  the  largest  number  of  stations 
given  event  has  been  i.  Usually  the  number 
and  data  irom  some  of  the  stations  have 


-  42  - 


never  been  seen  at  all.  The  situation  is  improving  through  the  efforts  of 

the  Albuquerque  Seismologtcal  Center,  which  now  has  responsibility  for 

operation  of  most  of  the  network,  but  their  predietion  is  tha  .  a  good 

quality  operation  is  still  some  time  in  the  future.  Instrument  malfunctions 

are  rather  frequent,  and  this  further  limits  the  comprehensiveness  of  the 
rirwi  lysis. 

Tin-  first  few  figure*  are  merely  examples  of  data  from  specific 
events.  They  will  show,  still  using  ns  a  descriptor,  what  can  be 
expected  for  various  sized  events,  both  explosions  and  earthquakes,  at 
various  distances.  Table  I  lists  parameters  of  the  events  shown,  together 
with  pertinent  station  data.  In  most  eases,  the  horizontal  components 
have  been  roughly  rotated  into  radial  and  transverse  components.  The 
transformation  is  necessarily  rough  because  the  frequency  responses 
and  nominal  magnifications  of  the  NS  and  EV.  components  at  each  station 
typically  clittci  cut,  and  in  Inc  absent  e  of  recent  data  the  reliability 
of  the  data  that  are  available  is  not  known.  The  instruments  were  treated 
as  if  the  horizontal  components  had  the  same  magnification  and  the  same 
frequency  response.  Except  in  one  case  the  results  are  quite  adequate, 
in  terms  of  separation  of  Love  and  Rayleigh  waves,  and  of  SH  and  SV 
type  'notion;  a  detailed  normalization  procedure  using  recent  calibration 
data  would  naturally  he  preferable.  Unfortunately  the  long-period  noise 
level  at  all  stations  in  the  network  is  much  higher  on  the  horizontal 
inst rumenl  s  than  „„  ,l,o  vertical  instruments.  Evidently  this  has  do 
piim.uily  with  tills,  as  the  vortical  instruments  are  relatively  unaffected. 
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tab 

L  E  I 

Date 

Origin  time 

Lat . 

1.011. 

Depth,  km 

rnb/MS 

Explosions  (| 

iresumed) 

Nov.  2. 

'72 

01  26  57.  6 

49.  9N 

78.  SE 

"0" 

6.  2  E 

Dec.  10, 

'72 

0-1  26  57.7 

49.  8N 

78.  OE 

"0" 

5.  7 

t  J 

"  27  08.  4 

50.  IN 

78.  8E 

"0" 

0-  0  (BRI<) 

E.  Kazakhstan 


Ea  rthguake  s 


'72 

14  03  16.  3 

52.  3N 

95.  4E 

"N"  5.5/4.  9 

RSFSR: 

Sayan  Mt.  Region 

'72 

00  14  10.  0 

35.  9N 

73.  4E 

55  5.  5 

Tadzhik  SSR 

'72 

06  59  59.  8 

36.  4N 

70.  6E 

204  4.  9 

Tadzhik  SSR 

and 

^  assigned  depth 

s  of  0  and 

j3  km.  r e s ne ctivf»  1  tj- 

Statistics  from  Preliminary  Detprmina*m„  ~c  ■ 

- lcrm*ndtion  of  Epicenters.  National 

Earthquake  Information  Service. 
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Figures  4.  ]  and  4.2  are  plots  of  seismograms  for  various 
presumed  explosions  and  earthquakes,  with  phases  identified  by  time 
of  arrival  and  polarization.  These  give  an  idea  of  the  signal-to-noise 
ratios  to  be  expected  for  various  event  sizes  and  types,  and  some  idea 
of  the  character  of  the  background  noise  on  the  records. 


Figures  4.  3  and  4.  4  concentrate  upon  the  excitation  of  P  and  PP 
waves  for  earthquakes  as  contrasted  with  explosions,  and  the  same  for 
vertically  and  horizontally  polarized  waves  of  the  S-type.  The  principal 
substantive  conclusion  from  this  study  up  to  this  time  is  that  then  only 
useful  long-period  body  waves  to  be  seen  from  the  presumed  explosions 
are  waves  of  the  SV  type,  ~.id  that  they  probably  originated  as  S  waves 
at  or  near  the  source.  This  subject  will  be  dealt  with  below. 

Figure  4.  5  illustrates  S-  waves  from  presumed  explosions  as 
recorded  on  several  instruments.  Even  though  the  signal-to-noise 
ratios  are  rather  low,  a  strong  case  can  be  made  for  the  appearance 
of  identical  wave  forms  on  the  various  records.  The  wave  form  as¬ 
sociated  with  SV-type  waves  can  be  seen  to  have  prograde  elliptical 
particle  motion  in  all  cases,  with  normal  dispersion.  Such  waves  fit  the 
definition  of  the  shear-coupled  PL  waves  described  by  Oliver  and  Major  (  ) 
and  as  in  the  case  of  earthquakes,  these  waves  are  more  prominent  on  the 
records  than  the  true  S  phase  itself.  It  is  noteworthy  that  a  phase  arriving 
at  the  time  of  PS  is  prominent  particularly  at  distances  beyond  90";  as  S 
is  lost  in  the  shadow  of  the  earth's  core.  PS  lakes  the  place  of  S  as  the 
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princip.il  body  phase,  again  showing  the  PL  wave  character.  If,  as 
seems  likely,  the  phase  is  actually  SP  (initial  long-period  P  waves 
being  unlikely  as  explained  previously),  then  the  energy  would  appear 
to  be  transmitted  primarily  as  P-  coupled  PL.  /dually  this  question 
is  indeterminate  at  present,  since  only  circumstantial  arguments  about 
the  generation  of  S  waves  are  available  on  the  basis  of  data  in  hand. 

Tiie  last  part  of  this  section  is  dedicated  to  a  discussion  of  a 
novel  theory  of  shear-wave  energy  generation  by  underground  explosions. 

It  is  an  ad-hoc  theory  developed  as  a  consequence  of  the  observed  long- 
period  SV  energy  from  presumed  explosions.  Many  of  tiie  predictions  of 
the  theory  are  similar  lo  those  of  the  fheory  that  S-  ”'avcs  arise  from 
conversion  of  P  waves  at  inhomogeneities  in  the  vicinity  of  the  explosion. 

A  test  of  the  theory  is  available,  however,  and  this  will  be  a  task  of  the 
immediate  future.  The  substance  of  the  theoretical  argument,  a  qualitative 
one  only  at  this  point,  follows. 

The  argument  as  originally  posed  relates  to  the  fact  that  the 
'lithostalic1  stress  state  in  a  rock  medium  under  gravitational  stress 
is  not  necessarily  a  hydrostatic  stress  state.  An  easy  way  to  see  this 
is  to  imagine  a  small  cylinder  of  rock.  When  the  cylinder  is  isolated, 
or  free  on  all  sides,  it  may  be  imagined  to  be  in  a  hydrostatic  stress 
state  at  zero  pressure,  neglecting  gravity  effects  within  the  cylinder. 
Locked -in  stresses  arc  not  considered.  L'ext,  consider  this  cylinder 
to  be  at  the  surface  of  a  half  space.  If  a  uniform  overburden  is  now 
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imposed  upon  the  half  space,  the  imaginary  cylinder  will  shorten  in  tht 
vertical  direction,  but  it  will  ho  prevented  from  freely  expanding  in  the 
lateral  direction  by  the  surrounding  ni.iicri.il.  It  is  then  clear  that  the 
cylinder  is  subjected  to  a  lateral  stress  at  depth  which  exceeds  the 
overburden  (vertical)  stress  component  by  the  factor  v  (1- v),  where  v 
is  Poisson's  ratio  for  the  material.  If  one  then  subtracts  the  hydrostatic 
pressure  cgh,  the  doviatoric  shear  stress  is  of  the  form  shown  in  figure 
4.  oa  that  is,  a  shearing  '-tress  of  the  order  of  30--40wo  of  the  overburden 
stress  is  likely  to  be  present  at  depth  in  the  earth  even  in  the  absence  of 
tile  so-called  "teconic  Stress"  which  has  been  invoked  to  explain  such 

2,  3 

things  as  the  Love  wave  radiation  observed  from  explosions  in  Nevada  ’  . 

This  is  not  to  say  that  that  type  of  tectonic  stress  may  not  be  present, 
b  it  rather  that  it  need  not  be  invoked  on  an  ad  hoc  basis  to  explain  observed 
shear  energy.  When  the  explosion  is  fired  it  immediately  forms  a  gas -filled 
cavity  several  meters  to  perhaps  a  hundred  meters  or  more  in  diameter. 
Naturally  the  gas  cannot  support  the  pre-existing  shear  stress  and  the 
surrounding  medium  must  adjust  elastically  in  shear  to  attain  equilibrium 
again.  Actually  the  adjustment  of  stress  amounts  to  a  collapse  in  shear  of 
the  cavity,  but  since  the  'cavity  collapse'  has  conic  to  denote  a  quite  dif¬ 
ferent  phenomenon  in  the  case  of  underground  explosions,  one  which  occurs 
at  least  minutes  or  hours  after  explosions,  this  deformation  will  be  termed 
''cavity  readjustment".  Ar.  originally  spherical  explosion  cavity  will 
be  transformed  to  an  oblate  sphe r ioidal  cavity  with  its  short  axis 
vertical,  consonant  with  the  lithoslatic  pre-stress.  See 
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figure  4.. ,b.Thv  . .  vvo. . .  , 

in  1  U(l  1  ur  lust.  .  t .< necmsly. 

The  ei;i  Stic  readjust  meat  of  the  neuly  ere  , ted  cntv  t 

)  vic.iuq  ij\ity  1 1 n  several 

interesting  Matures.  One  feature  i,  tha,  «.,  irs,  .notion  of  P  wave»  from 
.his  source  in  roughly  vertical  directions  nil,  bc  dU„a, in  contrast 

‘°  ,h°‘C  ^  lhC  eXpl0,i0n-  ?eC°"dl>’’  «“  «<■'«  -gion  for  the  seismic 
radiation  from  the  readjustment,  being  the  ton,  of  -hole  rock  surrounding 

•he  covi,y.  has  effective  radius  a,  leas,  as  large  as  tha,  of  the  equivalent  ’ 

cavity  due  to  the  espies, on  itself.  This  would  he  expected  affect  the 

frequency  spectrum  of  the  resulting  radiation.  ZnOnstie  offer,.,  including 

p  ast.e  o,  vis coela site  effects,  and  diffusion  of  g.,ses  through  the  fractured 

medium  surrounding  the  explosion  m.gh,  drag  ou,  the  completion  of  the 

Process,  which  would  undoubtedly  have  an  adverse  effect  upon  the  radiation 
of  higher  frequencies. 

Figure  4.  7  ind.cates  the  'focal  mechanism'  of  the  eiastic  readjustment 
expressed  in  terms  o,  firs,  motions.  At,  mentioned  above,  the  first  motion 

"fP  he  dilatations].  The  poiariration 

n,,<l  f'rSt  m"'i0"  "»  ^  ™  »c  predicted,  from  the  figure, 

for  , close i smic  arrival.,  and  the  surface  reflection  sS  would  arrive  a,  a 

. . *  ""■«  —  Civet,  that  the  symmetry  cf  the 

simple  problem  should  be  cvlindri,  il  ,1  v 

c>  Jindii.nl,  the  b  motion  should  bc  of  pure  SV 

polarisation.  Altliougl  in  vertic  il 

‘  ‘s  sccn°J1  tlle  mechanism  is  of  the 

double  couple  tvne  if  ic  n  r; 

-  •  figure  o,  revolution  about  the  vertical  axis,  wluelt 
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FIGURE  6  b 
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shows  the  collapse  or  the  crleinal H  coni>oncnt-  hewer  pal 

produced  instantaneously  hy  the  oxrlos?™11-"1  EaS-fiUed  c*vi, 

of  the  shear  modulus  within  the  cavity  f  r'fP°nso  fail, 

surface  is  shown  schematically  it  e2'i/'C  CollaPse  of  the  fi 

3t  should  he  expected. 

U-  First  ^tion  diagram  for  P  ,n.J  q  u 

r_e«id  justnieut  phenomenon  Cro-s*s  rhas*;s  °f  thc  cnvitv 

■about  the  vertical.  ’  action  of  a  figure  of  revolt, ti 


-  54  - 


Fig 


FREE  SURFACE 


4.7.  Typical  ray  paths  of  sS  and  S  phases  produced  by  the  cavity 

readjustment.  The  salient  point  is  that  S  and  sS  are  in  phase 

at  most  telesuismic  recording  stations,  in  contrast  to  the  case 
of  pP  and  P  from  both  the  explosion  and  cavity  readjustment, 
which  are  out  of  phase  and  hence  tend  to  cancel  one  another  at 
long  periods.  The  spectrum  of  sS  +  S  should  have  modulation 

appropriate  to  the  sS  -  S  time  with  a  maximum  at  zero  frequency. 

S  +  sS  should  have  substantial  energy  at  long  periods  as  is  observed. 
The  rosette  shows  the  amplitude  distribution  for  S.  Since  the 
elastic  cavity  readjustment  process  has  cylindrical  symmetry, 
the  drawing  shows  a  cross-section  of  a  figure  of  revolution  about 
the  vertical  axis. 
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pre-diets  th.,t  the  first  motion  patient  «ouM  be  a  function  of  distance  and 

FoH.nni.ely  a  . .  (,ropo,  js 

available,  since  it  predicts  that  the  amplitude  spectrum  of  the  'S'  phase 
Mould  have  a  maximum  at  aero  frequency,  and  a  first  minimum  at  1/(2  AT  ), 
uhere  At^  is  the  sS-S  time  interval  appropriate  to  the  depth  of  the  explosion. 
A  detailed  study  of  this  hypothesis  is  a  primary  goal  of  future  efforts. 

It  appears  that  the  project  to  study  long  period  body  waves  of 
underground  explosions,  has  produced  some  rather  intriguing  results. 

In  particular  the  long  period  S-  wave  radiation  from  explosions,  manifest 
primarily  in  particle  motion  of  the  PL  type  (normally  dispersed,  prograde 
elliptical,  SV  polarized)  is  surprisingly  prominent.  The  meager  obser¬ 
vational  data  support  a  novel  theoretical  model  which  includes  the  elastic 
readjustment  in  shear  of  the  medium  surrounding  the  explosion  cavity, 
without  the  necessity  of  the  type  of  'tectonic'  pre-stress  that  has  been 
invoked  in  the  past.  Predictions  of  first  motion  direction  and  pattern, 
including  independence  of  azimuth,  can  be  tested  and  should  be  valuable 
additions  to  the  suite  of  techniques  for  discrimination  of  explosions.  On 
the  other  hand.  Jong  period  body  waves  can  be  more  easily  studied  for 
known  earthquakes  of  substantially  smaller  magnitude  n^,  since  they  seem 
to  he  more  strongly  excited  by  earthquakes.  A  relative  excitation  criterion 
re seinhling  the  well  known  X-Wn^  criterion  appears  feasible  for  long 
period  body  waves  vs  in  (for  example). 
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5.  Classification  of  Asian  Earthquakes 
R.  Ganse  and  G.  Lundquist 

Daring  the  six  months  covered  by  this  report,  a  seismogram 
digitizing  system  and  associated  computer  software  have  been  established 
with  which  to  obtain  magnitudes,  focal  mechanisms,  and  seismic  source 
spectra  from  earthquake  seismograms.  Some  thirty  shallow  Asian  events 
were  selected  for  study  in  four  areas  bounded  by  40°  -  4  5°  N  and  70°  -  80°E; 
-193  -  3  5°  N  and  91 3  -  96c  E,  56°  -  38:N  and  94.  5°  -  97  5°  ,  and  50°  -  56°N  and 
104°-  1  13°E. 

Earthquale  Magnitudes. 

Three  separate  magnitudes  are  computed  from  every  available 
record  of  each  event.  Surface  wave  magnitudes,  M  ,  are  computed  from 
long-period,  vertical  records  using  the  International  Association  of  Seis¬ 
mology  and  Physics-  of  the  Earth's  Inter  or  (LASPEI)  formula. 
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where  the  amplitude  and  period  are  those  of  the  largest  wave  in  the  period 
range  1  7  to  23  seconds.  Long  period  body  wave  magnitudes  m^  are 
computed  according  to 
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where  amplitude  is  measured  from  peak  to  baseline  on  the  first  half  cycle 
of  the  P-  wave  motion  on  t lie*  long-period  vertical  component  seismograms. 


Q  (h.  A)  is  Gutenberg  .»nd  Richter's  well-known  empirical  correction 
Tor  distance  and  depth.  Short  period  body  wave  magnitude  rr^,  is  defined 


inb  *  lo"iol  <  Q 


amplitudi 


using  the  maximum  peak-io-peak  amplitude  found  in  the  first  3  £  cycles 
of  P-  wave  motion  on  the  short  period  vertical  seismogram.  In  addition, 
n\  an(1  MB  flre  to  be  computed  for  the  main  phase  of  any  events  judged  to 
be  a  multiple  event.  1  <>  date,  magnitudes  have  been  computed  for  about 
one-fourth  of  the  events  selected,  and  the  results  are  not  inconsistent 
with  worldwide  mb  -  M ^  relationships.  However,  results  are  much  more 
consistent  than  othc”  published  magnitudes  for  the  same  events,  probably 
clue  to  personal  reading  of  seismograms,  and  application  of  a  broader 
data  base  for  each  event. 

Preliminary  b  values  have  been  prepared  for  the  areas  under  study 
using  magnitudes  published  by  NEIC.  However,  the  number  of  events  is 

too  small  to  make  such  determinations  statistically  significant  in  three 
of  the  four  areas. 

Body  Wave  Focal  Vice ha n i t  ms. 

Focal  mechanisms  are  determined  from  P-  wave  first  motions, 
and  S-  wave  polarization  angles.  Some  P-  wave  first  motions  have  been 
gathered  from  the  JSC  Bulletin,  but  primary  emphasis  is  placed  upon 
readings  gathered  specifically  for  this  study.  S-  wave  horizontal  panicle 
motion  diagrams  are  constructed  for  all  S  waves  digitized  for  spectral 


analysis.  Experience  .so  far  indicates  that  one  or  two  good  S-  ave 

polarisation  angles  can  be  critically  important  in  disti-guishin'  betv.ecn 

possible  P-  wave  focal  mechanism  solutions.  A  catalogue  cl  predicted 

P-  and  S-  wave  motions  for  all  possible  double -couple  sources  (provided 

b>  Mi.  \ m.  Dillinger  ol  the  USGS)  will  be  used  to  aid  in  this  part  ol' the 
st  udy . 


tentative  focal  mechanisms  based  on  a  combined  P  and  S  solution 
fur  one  event  and  P  solutions  for  three  other  events  in  the  area  bounded 
by  29°  -  35°  N  and  9  1 r  -  9 ba  E  are  of  consistent  strike-slip  character.  The 
solution  will  be  refined  using  further  S-  wave  data  and  Russian  seismograms 
as  available.  One  combined  P  and  S  wave  focal  mechanism  in  the  area 
bounded  by  -10  -  45  X  and  t  (V  -  SO  E  agrees  with  the  predominantly  thrust 
focal  mechanisms  with  roughly  north-south  compression  published  by 

Molnur.-Cl L,  aL,  for  that  area.  Consistent  published  solutions  could  not 
be  found  for  the  other  area. 


So  i  s  m  i  c  S o u  3"  c  e  S  p c  c t  r  a . 

Body-wave  source  spectra  are  computed  from  both  short  and  long 
penod  V  WSSN  seismograms.  Seismograms  from  the  USSR  network  have 
been  ordered  and  spectra  will  also  be  derived  from  these.  Raw  P-  wave 


spectra  are  obtained  by  digitizing  ami  Fourie r -t ransforming  vertical 
component  records.  The  S-  wave  spectra  are  computed  from  the  trnns- 


verse  component  of  the  S-  wave. 


as  obtained  from  the  X-S  and  E-'U 


components  by  an  axis  rotation. 
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Thi*  IMW  spoct rn  ,11V  a»nvdt'(l  for  instrument  response  do  get 
station  g round -mot  ion  spectra)  and  then  for  anelaslic  attenuation.  The 
re  suiting  quasi-source  spectra  will  be  averaged  over  eight  or  more 
stations  with  good  azimuthal  distribution  to  give  the  fin'd  1  earthquake 
spectrum.  Seismic  moment,  stress  drop  and  fault  length  will  be  estimated 

from  spectra  of  sufficient  quality. 

Although  P-  wave  source  spectra  are  fairly  insensitive  to  the 

choice  of  Q,  the  locations  of  the  S-  wave  corner  frequencies  are  highly 
dependent  upon  the  anclastic  attenuation  model  used.  Separate  mecned 

quasi-source  spectra  were  constructed  by  applying  the  attenuation  versus 

frequency  corrections  predicted  by  the  velocity  model  CIT  11  CS2  and  Q 
model  QM  (low  Q)Z  and  CIT  208  (high  Q)3  model.  On  the  basis  of  com¬ 
patibility  of  the  P-  and  S-  wave  spectra,  results  from  four  events  seem 
to  favor  the  CIT  208  model.  An  example  of  the  spectra  for  the  different 
models  is  shown  in  Fie.  5.  1. 
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Mean  "Quasi-source  spectra"  (station  spectra  corrected  for 
anelastic  attenuation)  for  3  .April  1971  event.  Heavy  line  is 
spectrum  determined  from  long  period  seismograms;  light 
line  is  spectrum  determined  from  short  period  M'ismogrants. 
Diagram  A  is  P-  wave  spectrum  corrected  by  CIT  208  at¬ 
tenuation  model;  R  is  S-  wave  spectrum  corrected  by  CIT  208 
C  is  P-  wave  spectrum  corrected  by  CIT  1  1  CS2-QM;  D  is 
S-  wave  spectrum  corrected  by  CIT  1  1  CS2-QM. 


